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Abstract 

A search for a Higgs boson produced via vector-boson fusion and decaying into invisible 
particles is presented, using 20.3 fb”* of proton-proton collision data at a centre-of-mass 
energy of 8 TeV recorded by the ATLAS detector at the LHC. For a Higgs boson with a 
mass of 125 GeV, assuming the Standard Model production cross section, an upper bound 
of 0.28 is set on the branching fraction of // —> invisible at 95% confidence level, where 
the expected upper limit is 0.31. The results are interpreted in models of Higgs-portal dark 
matter where the branching fraction limit is converted into upper bounds on the dark-matter- 
nucleon scattering cross section as a function of the dark-matter particle mass, and compared 
to results from the direct dark-matter detection experiments. 
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1 Introduction 


Astrophysical observations provide strong evidence for dark matter (see Ref. [1] and the references 
therein). Dark matter (DM) may be explained by the existence of weakly interacting massive particles 
(WIMP) [2, 3]. The observed Higgs boson with a mass of about 125 GeV [4, 5] might decay to dark 
matter or neutral long-lived massive particles [6-10], provided this decay is kinematically allowed. This 
is referred to as an invisible decay of the Higgs boson [11-18]. 

This paper presents a search for invisible decays of a Higgs boson produced via the vector-boson fusion 
(VBF) process. In the Standard Model (SM), the process H —> ZZ —> 4v is an invisible decay of the 
Higgs boson, but the branching fraction (BF) is 0.1% [19, 20], which is below the sensitivity of the 
search presented in this paper. In addition to the VBF Higgs boson signal itself, there is a contribution 
to Higgs boson production from the gluon fusion plus 2-jets (ggF-i-2-jets) process, which is smaller than 
the VBF signal in the phase space of interest in this search. The ggF-t2-jets contribution is treated as 
signal. The search is performed with a dataset corresponding to an integrated luminosity of 20.3 fb“^ of 
proton-proton collisions at = 8 TeV, recorded by the ATLAS detector at the LHC [21]. 

The signature of this process is two jets with a large separation in pseudorapidity ^ and large missing 
transverse momentum^ The VBF process, in its most extreme topology (high dijet invariant mass 

for example), offers strong rejection against the QCD-initiated (lT,Z)-i-jets (V-i-jets) backgrounds. The 
resulting selection has a significantly better signal-to-background ratio than selections targeting the ggF 
process. 

The CMS Collaboration obtained an upper bound of 58% on the branching fraction of invisible Higgs bo¬ 
son decays using a combination of the VBF and ZH production modes [22]. Weaker limits were obtained 
using the Z(—> tt)H + F"™* ** signature by both the ATLAS and CMS collaborations [22, 23], giving upper 
limits at 95% CL of 75% and 83% on the branching fraction of invisible Higgs boson decays, respect¬ 
ively. By combining the searches in the Z{-^ {i)H and Z(^ bb)H channels, CMS obtained an upper 
limit of 81% [22]. Using the associated production with a vector boson, VH, where the vector-boson de¬ 
cays to jets and the Higgs boson to invisible particles, ATLAS set a 95% CL upper bound of 78% on the 
branching fraction o^ H ^ invisible [24]. Other searches for large in association with one or more 
jets were reported in Refs. [25-28]. These searches are less sensitive to Higgs-mediated interactions than 
the search presented here, because they are primarily sensitive to the ggF process and have significantly 
larger backgrounds. 

Assuming that the couplings of the Higgs boson to SM particles correspond to the SM values, global fits 
to measurements of cross sections times branching fractions of different channels allow the extraction of 
a limit on the Higgs boson’s branching fraction to invisible particles. The 95% CL upper limits on this 
branching fraction set by ATLAS and CMS are 23% and 21% respectively [29, 30]. There is an important 
complementarity between direct searches for invisible decays of Higgs bosons and indirect constraints 
on the sum of invisible and undetected decays. A simultaneous excess would confirm a signal, while a 
non-zero branching fraction of H ^ invisible in fhe global fif, buf no excess in fhe searches for Higgs 


* ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector 
and the z-axis along the beam direction. The x-axis points from the IP to the centre of the LHC ring, and the ^-axis points 
upward. Cylindrical coordinates (r, are used in the transverse plane, (p being the azimuthal angle around the beam direction. 
The pseudorapidity is defined as - In tan 0/2, where 6 is the polar angle. 

^ The transverse momentum is the component of the momentum vector perpendicular to the beam axis. 
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boson decays to invisible particles, would point toward other undetected decays or model assumptions as 
the source of the global fit result. 

In the search presented in this paper, the events observed in data are consistent with the background 
estimations. An upper bound on the cross section times the branching fraction of the Higgs boson decays 
to invisible particles is computed using a maximum-likelihood fit to the data with the profile likelihood- 
rafio fesf sfafisfic [31, 32]. A constrainl on fhe branching fracfion alone is obfained assuming fhe SM VBF 
and ggF producfion cross secfions, accepfances and efficiencies, for invisible decays of a Higgs boson 
wifh a mass mu - \25 GeV. 

In fhe confexf of models where dark mailer couples fo fhe SM particles primarily fhrough fhe Higgs 
boson [33], limifs on fhe branching fracfion of invisible Higgs boson decays can be inferprefed in WIMP- 
nucleon inferacfion models [34] and compared fo experimenfs which search for dark-maffer parficles via 
Iheir direcl inferacfion wifh fhe malerial of a delecfor [35-43]. 

The paper is organized as follows. The ATLAS delecfor is briefly described in Seclion 2. The modelling 
of fhe signal and background is presenled in Seclion 3. The dalasel, Iriggers, evenl reconslruclion, and 
evenl seleclion are described in Seclion 4. The background esfimafions are presenled in Seclion 5. In 
Seclion 6 , fhe sysfemafic uncerlainlies are discussed. The resulfs are shown in Seclion 7, and model 
inlerprefalions are given in Seclion 8 . Finally, concluding remarks are presenled in Seclion 9. 


2 Detector 

ATLAS is a mulfipurpose defector wifh a forward-backward symmelric cylindrical geomelry, described 
in delail in Ref. [21]. 

Al small radii from fhe beamline, fhe inner delecfor, immersed in a 2 T magnefic field produced by a Ihin 
superconducling solenoid located direclly inside fhe calorimeler, is made up of fine-granularily pixel and 
microslrip silicon deleclors covering fhe range I 77 I < 2.5, and a gas-hlled slraw-lube Iransilion-radialion 
Iracker (TRT) in fhe range \q\ < 2. The TRT complemenls fhe silicon fracker af larger radii and also 
provides elecfron idenlificalion based on Iransifion radiafion. The eleclromagnelic (EM) calorimeter is 
a lead/liquid-argon sampling calorimeter wifh an accordion geomelry. The EM calorimeler is divided 
into a barrel seclion covering I 77 I < 1.475 and Iwo end-cap secfions covering 1.375 < \'q\ < 3.2. Eor 
\r]\ < 2.5 if is divided info fhree layers in deplh, which are finely segmenled in rj and cf). An addifional 
Ihin presampler layer, covering \t]\ < 1 . 8 , is used fo correcf for fiucluafions in energy losses belween fhe 
producfion vertex and fhe calorimeter. Hadronic calorimefry in fhe region \t]\ <1.7 uses steel absorbers, 
and scinlillafor files as fhe aclive medium. Eiquid-argon calorimefry wifh copper absorbers is used in 
fhe hadronic end-cap calorimelers, which cover fhe region 1.5 < \r}\ < 3.2. A forward calorimeler using 
copper or lungsfen absorbers wifh liquid argon compleles fhe calorimeter coverage up fo \r]\ - 4.9. The 
muon speclromeler (MS) measures fhe curvalure of muon Irajecfories wifh \r]\ < 2.1, using fhree sfalions 
of precision drifl lubes, wifh calhode slrip chambers in fhe innermosl layer for \t]\ > 2.0. The defleclion is 
provided by a loroidal magnefic field wifh an integral of approximately 3 Tm and 6 Tm in fhe cenlral and 
end-cap regions of fhe ATEAS delecfor, respecfively. The MS is also insfrumenled wifh dedicated frigger 
chambers, namely resislive-plale chambers in fhe barrel and Ihin-gap chambers in fhe end-cap, covering 
\J]\ < 2.4. 
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3 Simulation 


Simulated signal and background event samples are produced with Monte Carlo (MC) event generators, 
and passed through a Geant4[44] simulation of the ATLAS detector [21, 45], or a fast simulation based 
on a parameterization of the response to the electromagnetic and hadronic showers in the ATLAS calor¬ 
imeters [46] and a detailed simulation of other parts of the detector and the trigger system. The results 
based on the fast simulation are validated against fully simulated samples and the difference is found to 
be negligible. The simulated events are reconstructed with the same software as the data. Additional pp 
collisions in the same and nearby bunch-crossings (pileup) are included by merging diffractive and non- 
diffractive pp collisions simulated with Pythia-8.165 [47]. The multiplicity distribution of these pileup 
collisions is re-weighted to agree with the distribution in the collision data. 

Both the VBF and ggF signals are modelled using Powheg-Box [48-53] with CTIO parton distribution 
functions (PDFs) [54], and Pythia-8.165 simulating the parton shower, hadronization and underlying 
event^. The VBF and ggF Higgs boson production cross sections and their uncertainties are taken from 
Ref. [55]. The transverse momentum (px) distribution of the VBF-produced Higgs boson is re-weighted 
to reflect electroweak (EW) radiative corrections computed by HAWK-2.0 [56]. These EW corrections 
amount to 10-25% in the Higgs boson px range of 150-1000 GeV. The ggE contribution to the signal is 
re-weighted [57, 58] so that the px distribution of the Higgs boson in events with two or more associated 
jets matches that of the next-to-leading-order (NEO) ggE-i-2-jets calculation in Powheg-Box MiNEO [59], 
and the inclusive distributions in jets match that of the next-to-next-to-leading-order (NNLO) and next- 
to-next-to-leading-logarithm (NNEE) calculation in HRes-2.1 [60, 61]. The effects of finite quark masses 
are also included [53]. 

The VT(—> ^y)-i-jets and Z(—> (£)+]ets processes are generated using Sherpa-L4.5 [62] including leading- 
order (EG) matrix elements for up to five partons in the final state with CTIO PDEs and matching these 
matrix elements with the parton shower following the procedure in Ref. [63]. The W{-^ fyj-i-jets and 
Z(^ processes are divided into two components based on the number of electroweak vertices 

in the Eeynman diagrams. Diagrams which have only two electroweak vertices contain jets that are 
produced via the strong interaction, and are labelled “QCD” Z-f-jets or W-i-jets. Diagrams which have 
four electroweak vertices contain jets that are produced via the electroweak interaction, and are labelled 
“EW” Z-f-jets or W-tjets [64]. The MC predictions of the QCD components of VT-i-jets and Z-rjets are 
normalized to NNEO in EEWZ [65, 66], while the EW components are normalized to VBENEO [67], 
including the jet px and dijet invariant mass requirements. The interference between the QCD and EW 
components of Z-rjets and W-f-jets is evaluated with Sherpa-L4.5 to be 7.5-18.0% of the size of the EW 
contribution depending on the signal regions. To account for this interference effect, the EW contribution 
is corrected with the estimated size of the interference term. Eigure 1 shows Eeynman diagrams for the 
signal and example vector-boson backgrounds. There are additional small backgrounds from tt, single 
top, diboson and multijet production. The tt process is modelled using Powheg-Box, with Pythia-8.165 
modelling the parton shower, hadronization and underlying event. Single-top production samples are 
generated with MC@NEO [68] for the s- and Wt-channel [69], while AcerMC-v3.8 [70] is used for 
single-top production in the t-channel. A top-quark mass of 172.5 GeV is used consistently. The AUET2C 
(AUET2B) [7 1] set of optimized parameters for the underlying event description is used for tt (single-top) 
processes, with CTIO (CTEQ6E1) [72] PDEs. Diboson samples WW, WZ and ZZ (with leptonic decays) 
are normalized at NEO and generated using Herwig-6.5.20 [73] with CTIO PDEs, including the parton 

^ The invisible decay of the Higgs boson is simulated by forcing the Higgs boson (with nijj = 125 GeV) to decay via H —> 
ZZ* ^ 4v. 
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Z+jets 


Figure 1: Example Feynman diagrams for the VBF //(—»invisible) signal and the vector-boson backgrounds. 


shower and hadronization, and Jimmy [74] to model the underlying event, whereas the WW, WZ, and ZZ 
{iqq, vvqq) processes are generated together with EW IF-i-jets and Z-i-jets samples. Diboson WW, 'WZ 
and ZZ (—> liqq, vvqq) samples generated using Sherpa-1.4.5 with CTIO PDFs and normalized to NLO 
in QCD [75] are used as a cross-check. Multijet and y-i-jet samples are generated using Pythia-8.165 with 
CTIO PDFs. 

4 Event selection 

The data used in this analysis were recorded with an F™®® trigger during periods when all ATFAS sub¬ 
detectors were operating under nominal conditions. The trigger consists of three levels of selections. The 
first two levels, FI and F2, use as inputs coarse-spatial-granularity analog (FI) and digital (F2) sums 
of the measured energy. In the final level, calibrafed clusfers of cell energies in fhe calorimefer [76] are 
used. Af each level, an increasingly sfringenf fhreshold is applied. The mosf stringenf requiremenf is 
^miss gQ QgY Because of furfher correcfions made in fhe offline reconstrucfed and the resolu¬ 

tions of the FI and F2 calculations, this trigger is not fully efficient until the offline is greater than 
150 GeV. 

Jets are reconstructed from calibrated energy clusters[77, 78] using the anti-A:, algorithm [79] with radius 
parameter R - 0.4. Jets are corrected for pileup using the event-by-event jet-area subtraction method [80, 
81] and calibrated to particle level by a multiplicative jet energy scale factor [77, 78]. The selected jets are 
required to have pT > 20 GeV and I 77 I < 4.5. To discriminate against jets originating from minimum-bias 
interactions, selection criteria are applied to ensure that at least 50% of the jet’s summed scalar track pT^ 
for jets within I 77 I < 2.5, is associated with tracks originating from the primary vertex, which is taken to 
be the vertex with the highest summed Pj of associated tracks. Information about the tracks and clusters 
in the event is used to construct multivariate discriminators to veto events with f7-jets and hadronic r-jets. 
The requirements on these discriminators identify fj-jets with 80% efficiency (estimated using tt events) 
[82-84], one-track jets from hadronic r decays with 60% efficiency (measured with Z ^ tt events), and 
multiple-track jets from hadronic r decays with 55% efficiency [85]. 

Electron candidates are reconstructed from clusters of energy deposits in the electromagnetic calorimeter 
matched to tracks in the inner detector [ 86 ]. Muon candidates are reconstructed by requiring a match 
between a track in the inner detector and a track in the muon spectrometer [87]. 

The selection defines fhree orfhogonal signal regions (SR), SRI, SR2a and SR2b. They are disfinguished 
primarily by fhe selecfion requiremenfs on fhe invarianf mass mjj of fhe fwo highesf-pT jets and fheir 
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separation in pseudorapidity Arjjj as shown in Table 1. The SRI selection requires events to have two 


Table 1; Summary of the main kinematic requirements in the three signal regions. 


Requirement 

SRI 

SR2a 

SR2b 

heading Jet pj 

>75 GeV 

>120 GeV 

>120GeV 

heading Jet Charge Fraction 

N/A 

>10% 

>10% 

Second Jet pp 

>50 GeV 

>35 GeV 

>35 GeV 


>1 TeV 

0.5 < rujj < 1 TeV 

> 1 TeV 

rjji X rjj2 


<0 

\^ rijj \ 

>4.8 

>3 

3 < \ Arjjj \ < 4.8 

\^< Pii \ 

<2.5 

N/A 

Third Jet Veto px Threshold 

30 GeV 

1 

>1.6 for 7 i, >1 otherwise 

>0.5 

rmiss 

>150GeV 

>200 GeV 


jets: one with pj > 75 GeV and one with pj > 50 GeV. The is constructed as the negative vectorial 
sum of the transverse momenta of all calibrated objects (identified electrons, muons, photons, hadronic 
decays of r-leptons, and jets) and an additional term for transverse energy in the calorimeter not included 
in any of these objects [88]. Events must have > 150 GeV in order to suppress the background from 
multijet events. To further suppress the multijet background, the two leading jets are required to have an 
azimuthal opening angle \A(pjj\ < 2.5 radians and an azimuthal opening angle with respect to the E™®® 
of |A(;(>y£.miss| >1.6 radians for the leading jet and |A(;(>y£.miss| > 1 radian otherwise. In the VBF process, 
the forward jets tend to have large separations in pseudorapidity (Arjjj), with correspondingly large dijet 
masses, and little hadronic activity between the two jets. To focus on the VBF production, the leading jets 
are required to be well-separated in pseudorapidity lArjjjl > 4.8, and have an invariant mass rrijj > 1 TeV 
Events are rejected if any jet is identified as arising from the decay of a f^-quark or a r-lepton. The rejection 
of events with fj-quarks suppresses top-quark backgrounds. Similarly, rejection of events with a r-lepton 
suppresses the VT(— > Tv)-i-jets background. Further, events are vetoed if they contain any reconstructed 
leptons passing the transverse momentum thresholds p^ > 10 GeV for electrons, pj>5 GeV for muons, 
or pj > 20 GeV for r-leptons. Finally, events with a third jet having pj > 30 GeV and \t]\ < 4.5 

are rejected. The SR2 selections are motivated by a search for new phenomena in final states with an 

energetic jet and large missing transverse momentum [25], and differ from those of SRI. First, the leading 
jet"^ is required to have pj > 120 GeV and \t]\ < 2.5. Additionally, the sub-leading jet is required to have 
Pt > 35 GeV, the Acpjj requirement is removed, the requirement on A^j^miss is relaxed to |A0y£.miss| > 0.5, 
and the E™'**® requirement is tightened to E™®® > 200 GeV. A common threshold of pj = 7 GeV is 
used to veto events with electrons and muons, and no r-lepton veto is applied. Finally in SR2, the E™®® 
computation excludes the muon contribution and treats hadronic taus like jets (this allows the modelling 
of IT-f-jets and Z-rjets in the control regions and signal regions using the same E™®® variable as discussed 
in Section 5). SR2 is further subdivided into SR2a with 500 < rrijj < 1000 GeV, pji x r]j 2 < 0, and 

\Ar]jj\ > 3, and SR2b with nijj > 1000 GeV, rjji x py 2 < 0 and 3 < \Ar]jj\ < 4.8. 


The “charge fraction" of this jet is defined as the ratio of the E/jt of tracks associated to the jet to the calibrated jet pi; this 
quantity must be at least 10% of the maximum fraction of the jet energy deposited in one calorimeter layer. The charged 
fraction requirement was shown to suppress fake jet backgrounds from beam-induced effects and cosmic-ray events [25]. 
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5 Background estimations 


In order to reduce the impact of theoretical and experimental uncertainties, the major backgrounds, Z ^ 
vv and W £v, are determined from measurements in a set of control samples consisting of Z ^ 
ox W ^ (v events = e/ju). In each of these control regions (CR), two additional jets are required, 
following the same requirements as the signal regions. The Z ^ control samples consist of events 
where the invariant mass of two same-flavour and opposite-sign leptons is consistent with the Z-boson 
mass, and so backgrounds in these control regions are small enough that they are taken from their MC 
predictions rather than from data-driven methods. In the W ^ iv control regions, the background from 
jets misidentified as leptons is more important, at least for the case of IT ^ ey. In SRI, the background 
from jets misidentified as leptons in the IT —> control regions is normalized using a fit that takes 

advantage of the distinctive shape of the transverse mass distribution 


/«T 


Ip^^E 


miss 

T 


1 - cos(A(;&^ ^miss ) 


( 1 ) 


of the lepton and and the charge asymmetry in IW~ events. In SR2, the background from jets 
misidentified as leptons in the IT —> /’v control regions is reduced by the requirements on nij and as 
discussed in Section 5.2. 


In order to use the control regions rather than the MC predictions for setting the IT-tjets and Z-rjets 
background normalizations, the MC predictions in each of the three signal regions and six corresponding 
Z(—> ee!pii)+]ots and 1T(—> ev//iv)-i-jets control regions are scaled by free parameters ki. There is one 
ki for each signal region and the corresponding control regions. In SRI for example, omitting factors 
that model systematic uncertainties, the expected number of events for Z{-^ vv)-tjets in the signal region 
is ZsRi = kiZ^, for Z{-^ i€)+]ot^ in the Z ^ control region is Zcr = kiZ^, and for W ^ tv 
in the VT —> control region is ITcr - kilT^. The scale factors ki are common for the Z-i-jets and 
IT-i-jets background normalizations. The scale factors ki are determined from the maximum Likelihood 
fit described in Section 7. The Z(—» ^fj-f-jets and the IT(^ ^y)-i-jets MC predictions thus affect the final 
estimates of Z(^ vv)-i-jets and IT(—> f’vj-rjets in the signal region through an implicit dependence on the 
MC ratio Zsr/Zcr and Zsr/VTcr for Z(^ yy)-f-jets, and VTsr/ITcr for 1T(^ f’yj-f-jets: 


ZsR ~ (Zsr/Zcr)mc ^ 

ZsR ~ (Zsr/M'cr)mc X W ^ f , (2) 

W^SR ~ (W^sr/W^cr)mc X 

Unique estimates of the Z(—> yy)-i-jets and 1T(—> f’yj-i-jets backgrounds in the signal region result from 
the simultaneous maximum likelihood fit to the control regions and signal region. 

The multijet background is estimated from data-driven methods as presented in Section 5.1. The data- 
driven normalizations for the Z-i-jets and IT-tjets backgrounds are described in Section 5.2. The back¬ 
ground estimations are validated in control regions with no signal contamination, and are in good agree¬ 
ment with observations in the validation control regions, as discussed in Section 5.3. In SRI and SR2, 
the smaller backgrounds of ti, single top and dibosons are taken from their MC predictions. 

Background contributions from the visible Higgs boson decay channels are suppressed by the signal 
region requirements described in Section 4. 
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5.1 Data-driven estimation of the multijet background 

Multijet events which have no prompt (from the primary interactions) neutrinos can pass the selec¬ 
tion due to instrumental effects such as the mis-measurement of the jet energy. Because of the very large 
rejection from the requirement, it is not practical to simulate this background, so it is estimated 
using data-driven methods instead. 

In the SR2 selections, the multijet background is estimated from data, using a jet smearing method as 
described in Ref. [89], which relies on the assumption that the of multijet events is dominated by 
fluctuations in the detector response to jets measured in the data. The estimated multijet background in 
SR2 is 24 ± 24 events (a 100% uncertainty is assigned to the estimate). 

In SRI, the multijet background is estimated from data as follows. A control region is defined where 
the £.miss requirement is inverted, so that the vector is in the direction of a jet in the event. The 
resulting sample is dominated by multijet events. The signal region requirements on the leading and 
sub-leading jet pj and on the trigger are applied as described in Section 4. The efficiency of each 
subsequent requirement is determined using this sample and assumed to apply to the signal region with 
the nominal A^^^-miss requirement. A systematic uncertainty is assessed based on the accuracy of this 
assumption in a control region with |A? 7 yy| < 3.8 and in a control region with three jets. To account for 
the A^ygmiss requirement itself, the A^yy requirement is inverted, requiring back-to-back jets in This 
sample is also multijet-dominated. Combining all the efficiencies with the observed control region yield 
gives an estimate of 2 ± 2 events for the multijet background in SRI. 


5.2 Estimations of the vv)-i-jets and ^v)-i-jets backgrounds 

To estimate the Z(^ vv)-i-jets background, both the Z(—» ee//r//)-f-jets and Vk(—> ev//iy)-i-jets control 
regions are employed. The IT-rjets background is estimated using IT(—> ev!pv)+]e.ts control regions. In 
the Z(—> ee)+iets control samples for SRI and IT(—> f’yj-i-jets control samples for SRI or SR2, electrons 
and muons are required to be isolated. Electron isolation is not required in the SR2 Z(—> ee)+]&ts control 
sample. For electrons, the normalized calorimeter isolation transverse energy, i.e. the ratio of the isolation 
transverse energy to lepton pT> is required to be less than 0.28 (0.05) for SRI (SR2), and the normalized 
track isolation is required to be less than 0.1 (0.05) within a cone AT? = - 0-3 for SRI 

(SR2). In the SRI selections, muons must have a normalized calorimeter isolation less than 0.3 (or < 0.18 
if pt < 25 GeV) and a normalized track isolation less than 0.12 within A/? = 0.3, whereas in the SR2 
selections, the scalar sum of the transverse momentum of tracks in a cone with radius 0.2 around the muon 
candidate is required to be less than 1.8 GeV. Electrons and muons are also required to point back to the 
primary vertex. The transverse impact parameter significance must be less than 3cr for both the electrons 
and muons, while the longitudinal impact parameter must be < 0.4(1.0) mm for electrons (muons). 

The Z(—> ee!pp)+]&t?, control regions are defined by selecting events containing two same-flavour, op¬ 
positely charged leptons with pj > 20 GeV and |m« - mz\ < 25 GeV, where mu and mz are the dilepton 
invariant mass and the Z-boson mass, respectively. In the control sample corresponding to the SRI se¬ 
lection, the leading lepton is required to have pj > 30 GeV. Triggers requiring a single electron or muon 
with pt > 24 GeV are used to select the control samples in SRI; in SR2, either a single-electron or 
trigger is used. The inefficiency of the triggers with respect to the offline requirements is negligible. In or¬ 
der to emulate the effect of the offline missing transverse momentum selection used in the signal region, 
the quantity is corrected by vectorially adding the electron (SRI and SR2) and muon transverse 



momenta (SRI only). All the Z(—> ee/ju/i)+jets events are then required to pass the other signal region 
selections. Backgrounds from processes other than Z(^ ee/ju/i)+jets are small in these control regions; 
the contributions from non-Z backgrounds are estimated from MC simulation. For Z —> ee (Z —> ////), the 
non-Z background is at a level of 1.6% (0.9%) of the sample. There is 50% uncertainty (mainly due to 
the limited numbers of MC events) on the non-Z background contamination in the Z control regions. The 
observed yield in the SRI Z control region, shown in Table 2, is larger than the expected yield by 16% 
but is compatible within the combined statistical uncertainties of MC simulation and data. In the SR2 

Table 2; Expected and observed yields for the SRI Z(—» ee/pyu)4-jets control sample in 20.3 fb ' of 2012 data. 
Expected contributions are evaluated using MC simulation, and the uncertainties are statistical only. 


Background 

SRI Z Confrol Regions 

Z(-^ ee)-i-jels Z{-^ pp)-i-jels 

QCD 

10.4 + 1.5 

14.0+ 1.5 


7.4 ± 0.8 

8.2 ± 0.8 

Of her Backgrounds 

0.3 ± 0.2 

0.2 ±0.1 

Tofal 

18.1 ± 1.7 

22.4 ± 1.7 

Dafa 

22 

25 


Table 3: Expected and observed yields for the SR2 Z(—> ee//ryu)4-jets control sample in 20.3 fb ' of 2012 data. 
Expected contributions are evaluated using MC simulation, and the uncertainties are statistical only. 


SR2 Z Confrol Regions 

SR2a 

SR2b 

Background 

Z(—> ec)+jels Z(^ pp)+jels 

Z(^ ce)+jels Z(^ pp)+jefs 

QCD z^e£ 

116 + 3 

121 ±4 

26 ± 2 

28 ±2 

Ewz^te 

17 ± 1 

17 ± 1 

16 ± 1 

16 ±2 

Ofher backgrounds 

8 ± 1 

10 ±2 

2 ± 1 

3 ± 1 

Tofal 

141 ± 3 

148 ±5 

44 ± 3 47 ± 3 

Dafa 

159 

139 

33 

38 


control regions, the observed and expected yields differ by 10% as shown in Table 3 but are compatible 
within the total statistical and systematic uncertainties (see Section 6). The emulated distributions 
for the Z control regions are shown in Figures 2 and 3 for SRI and SR2 respectively. Because the muon 
momentum is excluded from the definition in SR2, the "emulated" label is omitted from Figure 3(b). 

The VF(^ ev/yuy)-i-jets control regions are similarly defined by selecfing evenfs confaining one lepfon 
wifh fransverse momenfum px > 30 GeV (25 GeV) in fhe case of SRI (SR2), and no additional leptons 
wifh pt > 20 GeV. The is emulafed in fhe same way as for fhe Z —> cc/pp confrol region and 
evenfs are required fo pass fhe signal region selections on jefs and In SRI, fhe confribufions of fhe 
fhree lepfon flavours fo fhe fofal W ^ (v background affer all fhe requiremenfs are 20% for W —> ey, 
20% for IT ^ py, and 60% for W —> tv. The equal confribufions from W ^ ev and IT ^ py suggesf 
fhaf fhese are evenfs where fhe lepfon is below ifs px fhreshold or suflicienfly far forward fo escape fhe 
jef veto, and nol evenfs where fhe lepfon is misidenfified as a fag jef (since muons deposif liffle energy 
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in the calorimeter and would therefore not be identified as a jet). This expectation is checked explicitly 
for the case of VT —> tv, by using MC truth information about the direction of the r-lepton to find the A/? 
between the r-lepton and the nearest reconstructed tag jet. The component with A/?/,- < 0.4 is completely 
negligible after the signal region requirements, indicating that the lepton tends to be recoiling against 
the tag jets rather than being aligned with them. For SRI, four W control regions are considered using 
different charge samples for W‘/VT“ ^ evliuv since 1T(—> fvj-i-jets is not charge symmetric as shown in 
Table 4, whereas in SR2, only two control regions 1T(—> ey//rv)-i-jets are used as shown in Table 5. 


Table 4: Expected and observed yields for the SRI W iv control sample, after all requirements in 20.3 fb ' 
of 2012 data. The multijet background is estimated using the data-driven method described in the text; all other 
contributions are evaluated using MC simulation. Only the statistical uncertainties are shown. 


Background 

IT+ ^ ev 

SRI W Conlrol Regions 
IT- ^ ev IT+ ^ pv 

IT- 

pv 

QCD W ^{v 

92.3 ± 7.2 

55.1 ±5.3 

85.5 ± 7.0 

43.8 

±4.6 

EWW ^{v 

99.4 ± 4.0 

52.5 ± 2.9 

81.9 ± 3.7 

39.1 

±2.5 

QCD 

3.4 ± 0.6 

4.4 ± 0.9 

6.4 ± 1.1 

5.0 

±0.9 


2.5 ± 0.3 

2.9 ± 0.3 

2.7 ± 0.3 

3.2 

±0.3 

Mullijel 

28.0 + 6.8 

28.0 ± 6.8 

1.6 ±2.6 

1.6 

±2.6 

Olher backgrounds 

4.0 ± 0.7 

1.8 ±0.4 

3.2 ± 0.7 

1.0 

±0.3 

Tolal 

230 ± 11 

145 ±9 

181 ± 8 

93.7 ± 5.9 

Dala 

225 

141 

182 

98 


Table 5: Expected and observed yields for the SR2 VE(—» ev/pv)4-jets control sample in 20.3 fb ' of 2012 data. Ex¬ 
pected contributions are evaluated using MC simulation, and the uncertainties are statistical only. The discrepancy 
in the Vk(—> pv)H-jets SR2b control region is due to a mis-modelling of the W pj. The agreement improves when 
the systematic uncertainties (discussed in Section 6) are included. 

SR2 W Control Regions SR2a SR2b 


Background 

1T(^ ey)±jels 

IT(^ juy)±jels 

IT(^ ey)±jels 

1T(^ yuy)±jels 

QCD W 

595 ± 12 

906 ± 15 

122 ±5 

201 ±7 

EWW ^£v 

149 ±5 

214 ± 6 

121 ±4 

184 ±5 

QCD z^a 

5.8 ±0.9 

23 ± 1.6 

1.6 ± 0.4 

4.5 ± 0.6 

EWZ^e£ 

0.4 ±0.1 

0.5 ± 0.2 

2.0 ± 0.4 

2.7 ± 0.8 

Mullijel 

13 ±3 

0±0 

3 ± 1 

0±0 

Olher backgrounds 

44 ±4 

78 ±7 

13 ±2 

19 ±3 

Tolal 

807 ± 14 

1222 ± 18 

263 ±7 

411 ±9 

Dala 

783 

1209 

224 

295 


In the W{-^ ev/pv)+iets control regions corresponding to the SRI selection, a fit to the transverse mass 
defined in Eq. (1) is used lo eslimale Ihe mullijel background. In order lo oblain an explicif measuremenl 
and uncerfainfy for Ihe background from mulfijels, no requiremenls are made on E™®® and mv. Because 
Ihe mullijel background does no! have a prompl neufrino, Ihe lends lo be lower and lo poinl in Ihe 


11 



















direction of the jet that was misidentified as a lepton. As a result, the multijet background tends to have 
significantly lower mj than the IV+jets contribution. Control samples modelling the jets misidentified as 
leptons in mulfijef evenfs are consfrucfed by selecting evenfs fhaf pass fhe ik+jefs confrol region selecfion, 
excepf for cerfain lepfon idenfificafion criferia: for elections, some of fhe EM calorimefer shower shape 
requiremenfs are loosened and fully identified elections are removed, while for muons, fhe fransverse im- 
pacf parameter (do) requiremenf which suppresses muons originafing from heavy-flavour jefs is reversed. 
To obfain fhe normalizafion of fhe mulfijef background in fhe IT-i-jefs confrol region, femplafes of fhe 
niT disfribufion for processes wifh prompf lepfons are faken from MC simulafion. Shape femplafes for 
fhe backgrounds from mulfijef evenfs are consfrucfed by summing fhe observed yields in confrol samples 
obfained by inverting fhe lepton identification and do requiremenfs, and subfracfing fhe expected contii- 
bufions from IT-i-jefs and Z-i-jefs evenfs using MC. Since fhe misidenfified-jef samples are expected to be 
charge-symmefric, fhe same shape femplafe and normalizafion factor is used fo model bofh charge caf- 
egories of a given lepton flavour {e or /i). To determine fhe 1T(—> ^v)-i-jels background normalizafion, a 
fif fo fhe fransverse mass mj of fhe lepfon and is used. The VT(—> fvj-f-jefs confribufion, however, is 
nof charge symmefric, so fhe differenf charge samples are kepi separale in fhe simullaneous fil fo four mj 
disfribulions, one for each lepton flavour and charge combinafion shown in Figure 4. There are fhree free 
normalizations in fhe fif: one for evenfs wifh a prompf lepfon, one for evenfs where a jel is misidenfified 
as an election, and one for evenfs where a jel is misidenfified as a muon. The normalization factor for the 
prompt leptons in the /wt fit is 0.95 + 0.05 (stat). 

In the W —> ev control region corresponding to the SR2 selections, the background from multijet events 
is rejected by requiring that the (corrected by vectorially adding the election transverse momentum) 
be larger than 25 GeV and that the transverse mass be in the range 40 < mj < 100 GeV. The selected 
electron is required to pass both the track and calorimeter isolation requirements. The tight requirements 
on electron isolation and Ej'^^ greatly reduce the multijet background relative to the other backgrounds. 
The residual multijet background in the IT —> ev control region is at the level of 1% of the total control 
region background, with an uncertainty of 100%. For the W ^ fiv control region corresponding to 
the SR2 selections, the selected muon is required to pass only the track isolation requirement and the 
transverse mass is required to be in the range 30 < /wt < 100 GeV. An attempt is made to estimate the 
residual multijet background in the W ^ jiv control region using a control sample with inverted muon 
isolation. The residual background from multijet events is negligible. Figure 5 shows the mx distributions 
in the SR2 IT-rjets control regions. 


5.3 Validation of Data-Driven Estimations 

To validate the background estimates for SRI, two signal-depleted neighbouring regions are defined by 
(1) reversing fhe veto againsf fhree-jef evenfs and requiring fhaf fhe fhird jef in fhe evenf has fransverse 
momenfum > 40 GeV, and (2) reversing bofh fhe jef veto wifh a > 30 GeV requiremenf and fhe 
jef rapidify gap wifh a |A? 7 /y| < 3.8 requiremenf. Good agreemenf befween expecfafion and observafion is 
found in fhese validation regions, as shown in Table 6. 


6 Systematic uncertainties 


The experimenfal uncerfainfies on fhe MC predicfions for signals and backgrounds are dominated by 
uncerfainfies in fhe jef energy scale and resolufion [77]. This includes effecfs such as fhe p dependence 
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Figure 4: The transverse mass distributions used in the SRI VT+jets control regions after all requirements except 
for the > 150 GeV requirement: (a) e^v, (b) —* e^v, (c) and (d) W~ yu'v. 
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Figure 5: The transverse mass distributions in the SR2 VT-i-jets control regions after all requirements: (a) VT —> ev 
and (b) W ^v. 
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Table 6: Expected and observed yields for the validation regions in 20.3 fb ' of data. 3-jet: reversal of the veto 
against three-jet events by requiring > 40 GeV; and 3-jet and lAr/jyl < 3.8: requirements of both lAr/yjl < 3.8 
and pj > 30 GeV. Contributions from fT-rjets and Z-i-jets are normalized to data-driven estimates. The fT-tjets and 
Z-i-jets uncertainties include MC statistics from both the selected region and the corresponding control region, and 
the number of data events in the control regions. The other numbers are evaluated using MC simulation and their 
uncertainties indicate only statistical uncertainty. 


Process 

3-jef 

3-jef and lAp,/! < 3.8 

ggF signal 

6.2 ±3.1 

- 

VBF signal 

19.9 ± 1.4 

4.7 ±0.6 

Z{-^ yy)-i-jels 

97 ± 10 

111 ± 10 

VF(^ ^y)-i-jels 

78.5 ±6.5 

73 ± 10 

Mulfijef 

19.9 ±21.8 

- 

Ofher backgrounds 

2.2 ± 0.3 

0.5 ±0.1 

Tofal 

198 ± 25 

185 ± 14 

Dafa 

212 

195 


of the energy scale calibration and the dependence of the energy response on the jet flavour composition, 
where flavour refers to the gluon or light quark initiating the jet. Uncertainties related to the lepton 
identification in the control regions and lepton vetoes are negligible. Luminosity uncertainties [90] are 
applied to the signal and background yields that are obtained from MC simulation. 

Theoretical uncertainties on the IV-i-jets and Z-i-jets contributions to both the signal and control regions 
are assessed using Sherpa, and cross-checked with MCFM [75] and VBFNLO [67] for the EW and QCD 
processes respectively, and by a comparison between Sherpa and Alpgen [91] for the latter process. In 
all cases, the uncertainties are determined by independently varying the factorization and renormaliza¬ 
tion scales by factors of 2 and 1 /2, keeping their ratio within 0.5-2.0. The parton distribution function 
uncertainties are evaluated with the CTIO error sets [54]. The uncertainty on the ggF yield due to the 
jet selection is evaluated using Stewart-Tackmann method [92]. Uncertainty in the px distribution of the 
Higgs boson in ggF is evaluated from scale variations in HRes following the re-weighting of the px distri¬ 
bution [60, 61] as mentioned in Section 1. To assess the level of theoretical uncertainty on the jet veto, the 
variation in the predicted VBF cross section with respect to shifts in the renormalization and factorization 
scales as well as with respect to uncertainty in the parton-shower modelling is measured using Powheg- 
Box NLO generator matched to Pythia and to Herwig. The effect of the parton shower on the QCD 
IF-i-jets and Z-f-jets background estimations is obtained by comparing simulated samples with different 
parton shower models. As shown in Table 7, where the main systematic uncertainties are summarized, 
using the MC predictions of Zsr/IFcr and IFsr/W^cr ratios reduces the systematic uncertainties in the 
final Z-i-jels and VF-i-jels background esfimafes. The Z(—» ^Q-i-jels/VF(—> 7’y)-i-jels rafio is checked in dafa 
and MC, and no discrepancy larger fhan 10% is observed, consisfenf wifh fhe residual fheory uncerfainfies 
on fhe Zsr/VFcr ratios shown in Table 7. 
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Table 7: Detector and theory uncertainties (%) after all SR or CR selections. For each source of uncertainty, where 
relevant, the hrst and second rows correspond to the uncertainties in SRI and SR2 respectively. The ranges of 
uncertainties in the Z or IT column correspond to uncertainties in the Z+jets and IT+jets MC yields in the SR or 
CR. The search uses the uncertainties in the ratios of SR to CR yields shown in the last column. 


Uncertainty 

VBE 

ggF 

Zor IT 

Zsr/ITcr or ITsr/ITcr 

Jet energy scale 

16 

9 

43 

12 

17-33 

0-11 

3-5 

1 ^ 

Jet energy resolution 

Negligible 

3.1 

Negligible 

3.2 

Negligible 

0.2-7.6 

Negligible 

0.5-5 .8 

Euminosity 

2.8 

2.8 

2.8 

Irrelevant 

QCD scale 

0.2 

7.8 

5-36 

7.5-21 

7.8-12 

1-2 

PDE 

2.3 

2.8 

7.5 

3-5 

0 . 1 - 2.6 

1-2 

Parton shower 

4.4 

41 

29 

9-10 

5 

Veto on third jet 

Negligible 

Negligible 

Higgs boson pj 

Negligible 

9.7 

Irrelevant 

Irrelevant 

MC statistics 

2 

0.6 

46 

13 

2.3-6.4 

0.8^.5 

3.3-6.6 


7 Results 

Figures 6 and 7 show the and the mjj distributions after imposing the requirements of SRI and 
SR2 respectively. There is good agreement between the data and the background expectations from the 
SM, and no statistically significant excess is observed in data. The limit on the branching fraction of 
H invisible is computed using a maximum-likelihood fit to the yields in the signal regions and the iT(^ 
ev/fj.v)+iets and Z(—> ee////r)-i-jets control samples following the CL 5 modified frequenfisf formalism [93] 
wifh a profile likelihood-rafio fesf sfafisfic [31]. Expecfed signal and background disfribufions in fhe signal 
and confrol regions are defermined from MC predicfions, wifh fhe excepfion of fhe mulfijef backgrounds, 
which use fhe dafa-driven mefhods described in Secfion 5. Sysfemafic uncerfainfies are parameferized 
as Gaussian consfrained nuisance paramefers. The nuisance paramefer for each individual source of 
uncertainly is shared among fhe expecfed yields so lhal its correlated effect is taken into account. The 
relative weight of the Z(—> ee/ju/i)-i-jets and W(-^ ev/yuv)-tjets in the control regions is determined by the 
maximization of the likelihood function. 

One global likelihood function including all three signal regions and the six corresponding control regions 
is constructed with only the signal yields and correlated uncertainties coupling the search regions. The 
theoretical uncertainties are taken to be uncorrelated between the EW and QCD processes and uncorrel¬ 
ated with the scale uncertainty on the signal. The uncertainties which are treated as correlated between 
the regions are: 
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Figure 6; Data and MC distributions after all the requirements in SRI for (a) and (b) the dijet invariant mass 
nijj. The background histograms are normalized to the values in Table 8. The VBF signal (red histogram) is 
normalized to the SM VBF Higgs boson production cross section with BF{H —»invisible) = 100%. 



(a) distribution 


(b) irijj distribution 


Figure 7; Data and MC distributions after all the requirements in SR2 for (a) and (b) the dijet invariant mass 
nijj. The background histograms are normalized to the values in Table 8. The VBF signal is normalized to the SM 
VBF Higgs boson production cross section with BF(H —» invisible) = 100%. 
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• Uncertainty in the luminosity measurements. This impacts the predicted rates of the signals and the 
backgrounds that are estimated using MC simulation, namely ggF and VBF signals, and tt, single 
top, and diboson backgrounds. 

• Uncertainties in the absolute scale and resolution of the reconstructed jet energy. 

• Uncertainties in the modelling of the parton shower. 

• Uncertainties in renormalization and factorization scales. 

Table 8 shows signal, background and data events after the global fit including the effects of systematic 
uncertainties, MC statistical uncertainties in the control and signal regions, and the data statistical uncer¬ 
tainties in the control regions. The post-fit values of the Z-i-jets and IF-i-jets background normalization 
scale factors k,, discussed in Section 5, are 0.95 ± 0.21, 0.87 ± 0.17 and 0.74 ± 0.12 for SRI, SR2a and 
SR2b and their control regions, respectively. As shown in Table 8, the signal-to-background ratio is 0.53 
in SRI, and 0.09 and 0.21 in SR2a and SR2b respectively, for BF{H invisible) = 100%. Fits to the 

Table 8; Estimates of the expected yields and their total uncertainties for SRI and SR2 in 20.3 fb ' of 2012 data. 
The Z(—> vv)4-jets, W{^ f’v)4-jets, and multijet background estimates are data-driven. The other backgrounds 
and the ggF and VBF signals are determined from MC simulation. The expected signal yields are shown for 
ififf - 125 GeV and are normalized to BF(// —> invisible) = 100%. The fk-i-jets and Z-i-jets statistical uncertainties 
result from the number of MC events in each signal and corresponding control region, and from the number of data 
events in the control region. 


Signal region 

Process 

SRI 

SR2a 

SR2b 

ggF signal 

20±15 

58± 22 

19± 8 

VBF signal 

286±57 

182± 19 

105±15 

Z(^ yy)-l-jets 

339±37 

1580± 90 

335±23 

W{-^ ^y)-i-jets 

235±42 

1010± 50 

225±16 

Multijet 

2 ± 2 

20 ± 20 

4±4 

Other backgrounds 

1±0.4 

64± 9 

19±6 

Total background 

511 ±61 

2680±130 

583±34 

Data 

539 

2654 

636 


likelihood function are performed separately for each signal region and their combination, and the 95% 
CL limits on BF(// ^ invisible) are shown in Table 9. 


Table 9; Summary of limits on BF{H —> invisible) for 20.3 fb ' of 8 TeV data in the individual search regions and 
their combination, assuming the SM cross section for oth = 125 GeV. 


Results 

Expected 

-i-lo- 

-Icr 

-l-2cr 

-2cr 

Observed 

SRI 

0.35 

0.49 

0.25 

0.67 

0.19 

0.30 

SR2 

0.60 

0.85 

0.43 

1.18 

0.32 

0.83 

Combined Results 

0.31 

0.44 

0.23 

0.60 

0.17 

0.28 
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The agreement between the data and the background expectations in SRI is also expressed as a model- 
independent 95% CL upper limit on the fiducial cross section 

O'M = CT X BF X J?l, 

N 

where the acceptance J?1 is the fraction of events within the fiducal phase space defined af fhe MC frufh 
level using fhe SRI selections in Secfion 4, N fhe accepfed number of evenfs, X the integrated luminosity 
and e the selection efficiency defined as fhe rafio of selecfed evenfs fo fhose in fhe fiducial phase space. 
Only fhe sysfemafic uncerfainfies on fhe backgrounds and fhe infegrafed luminosify are faken info accounf 
in fhe upper limif on ctm, shown in Table 10. In SRI, fhe accepfance and fhe evenf selection efficiency, 


(3) 

(4) 


Table 10: Model-independent 95% CL upper limit on the fiducial cross section for non-SM processes 0-^4 in SRI. 


SRI 

Expected 

-l-lcr 

-Icr 

+2(r 

-2cr 

Observed 

Fiducial cross section [fb] 

4.78 

6.32 

3.51 

8.43 

2.53 

3.93 


estimated from simulated VBF H —> ZZ ^ 4v events, are (0.89 ± 0.04)% and (94 ± 15)% respectively. 
The uncertainties have been divided such that the theory uncertainties are assigned to the acceptance and 
the experiment uncertainties are assigned to the efficiency. 


8 Model interpretation 


In the Higgs-portal dark-matter scenario, a dark sector is coupled to the Standard Model via the Higgs 
boson [9, 10] by introducing a WIMP dark-matter singlet that only couples to the SM Higgs doublet. 
In this model, assuming that the dark-matter particle is lighter than half the Higgs boson mass, one would 
search for Higgs boson decays to undetected (invisible) dark-matter particles, e.g. H —> XX- The upper 
limits on the branching fraction to invisible particles directly determine the maximum allowed decay 
width to the invisible particles 


F 


inv 

H 


BF{H —> invisible) 

1 - BF{H invisible) 


(5) 


where F// is the SM decay width of the Higgs boson. Adopting the formulas from Ref. [10], the decay 
width of the Higgs boson to the invisible particles can be written as 


T-inv 
^ H^SS 


T-inv 
^ H^VV 


pinv 


6AnmH 
256mny \ 
327rA2 ’ 


-t 12 


4 

my 

4 


m 


H 


( 6 ) 

(V) 

( 8 ) 


for the scalar, vector and Majorana-fermion dark matter, respectively. The parameters Ahss^ Ahvv, 
Auffl^ are the corresponding coupling constants, v is the vacuum expectation value of the SM Higgs 
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doublet, = y 1 (x = S ,'V, f), and is the WIMP mass. In the Higgs-portal model, the 

Higgs boson is assumed to be the only mediator in the WIMP-nucleon scattering, and the WIMP-nucleon 
cross section can be written in a general spin-independent form. Inserting the couplings and masses for 
each spin scenario gives: 


cr 


SI 

SN 


cr 


si 

VN 


cr 


SI 

/W 


'^HSS "‘N-IN 

1 Gnmjj {ms + 

A^ 

'^HVV '"nJn 

1 {my + mjsf)^ 

^4 2 f2 

{mf + miy)^ 


(9) 

( 10 ) 

( 11 ) 


where is the nucleon mass, and is the form factor associated to the Higgs boson-nucleon coupling 
and computed using lattice QCD [10]. The numerical values for all the parameters in the equations above 
are given in Table 11 . 


Table 11: Parameters in the Higgs-portal dark-matter model. 


Vacuum expectation value 

r/V2 

174 GeV 

Higgs boson mass 

mn 

125 GeV 

Higgs boson width 

Th 

4.07 MeV 

Nucleon mass 

mN 

939 MeV 

Higgs-nucleon coupling form factor 

In 

0 33‘^°-^‘’ 

^•-^■^-0.07 


The inferred 90% CL branching fraction limit for H invisible, translated into an upper bound on 
the scattering cross section between nucleons and WIMP, is shown in Figure 8 compared to the results 
from direct detection experiments. The WIMP-nucleon cross-section limits resulting from searches for 
invisible Higgs boson decays extend from low WIMP mass to half the Higgs boson mass, and are com¬ 
plementary to the results provided by direct detection experiments that have limited sensitivity to WIMP 
with mass of the order of 10 GeV and lower [35, 37-41, 43]. This is expected as the LHC has no limit¬ 
ations for the production of low-mass particles, whereas the recoil energies produced in the interactions 
of sub-relativistic WIMP with nuclei in the apparatus of a direct detection experiment are often below 
the sensitivity threshold for small WIMP masses. The aforementioned correlation between the branching 
fraction of Higgs boson decays to invisible particles and the WIMP-nucleon cross section is presented in 
the effective field theory framework, assuming that the new physics scale is (9(afew) TeV, well above the 
scale probed at SM Higgs boson mass. Adding a renormalizable mechanism for generating the fermion 
and vector WIMP masses could modify the correlation between the WIMP-nucleon cross section and the 
branching fraction of Higgs boson decays to invisible particles [94]. 


9 Conclusions 

A search for Higgs boson decays to invisible particles is presented. The search uses data events with 
two forward jets and large missing transverse momentum, collected with the ATLAS detector from 20.3 
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Figure 8; The WIMP-nucleon cross section as a function of the WIMP mass. The exclusion limits [35-39] of the 
direct detection experiments are compared to the ATLAS results from the BF{H —» invisible) limit in the Higgs- 
portal scenario, translated into the WIMP-nucleon cross section using the formulas from Ref. [10]. The exclusion 
limits are shown at 90% CL. The error bands on the ATLAS results indicate the uncertainty coming from the 
different estimations of the Higgs-nucleon coupling form factor [95, 96]. 


fb“^ of pp collisions at = 8 TeV at the LHC. Assuming the SM production cross section, acceptance 
and efficiency for invisible decays of a Higgs boson with a mass of 125 GeV, a 95% CL upper bound is 
set on the BF(// ^ invisible) at 0.28. The results are interpreted in the Higgs-portal dark-matter model 
where the 90% CL limit on the BF(// ^ invisible) is converted into upper bounds on the dark-matter 
nucleon scattering cross section as a function of the dark-matter particle mass. The ATLAS limits are 
complementary to the results from the direct dark-matter detection experiments. 
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